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Abstract 

We use Pade Approximants to obtain improved predictions for the anomalous 
magnetic moments of the electron and the muon. These are needed because 
of the very precise experimental values presently obtained for the electron, 
and soon to be obtained at BNL for the muon. The Pade prediction for the 
QED contribution to the anomalous magnetic moment of the muon differs 
significantly from the naive perturbative prediction. 
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Two of the most important tests of quantum electrodynamics (QED) are 
the comparisons between theory and experiment of the anomalous magnetic 
moments of the electron and the muon, a e and a M respectively, where a = 
(g—2)/2. The latest Penning trap measurements of the electron and positron 
anomalies obtained by the University of Washington groupS are: 

a ex P t = 115 9 65 2188.4(4.3) x KT 12 (1) 

and 

a e X = 1159652187.9(4.3) x 1(T 12 (2) 

The figures in brackets represent the error in the last 2 figures, a convention 
we will follow throughout this paper. Taking the average of eqs (H) and @, 
one finds 

afP* = 1159652188.2(3.0) x 1(T 12 (3) 

The most accurate measurement for the muon anomaly comes from the 
CERN g— 2 experimental in which it was found that 

a^ x _ pt = 1165936(12) x 1(T 9 (4) 

and 

a; x + pt = 1165910(11) x 10~ 9 (5) 
and the combined result is 



a 



,)f pt = 1165923(9) x 10~ 9 (6) 



where correlations are taken into account in combining the errors. A new 
g—2 muon experiment is being done at Brookhaven National Laboratory 
(BNL)l, and an improvement in the accuracy by a factor of about 20 is 
expected. In order to compare properly theory and experiment, one must 
improve correspondingly the accuracy of the theoretical predictions. 

In an heroic feat, Kinoshitai has calculated a e in eighth order and Ki- 
noshita, Nizic, OkamotoH and MarciancS have calculated a M in eighth or- 
der. Moreover, there have been some recent improvements in the analytic 
calculations!^ of a e and a M . 

There have recently been several papers estimating coefficients in Per- 
turbative Quantum Field Theory (PQFT) using Pade ApproximantslEl. 
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This procedure is known to give significant improvements on naive perturba- 
tive calculations in many condensed-matter applications^], removes a large 
part of the discrepancy between experiment and QED calculations of the 
ortho-positronium decay rate0'0 and agrees well with other estimates of 
higher-order perturbative coefficients in QCdH'0. 

In this paper we will use Pade Approximants (PA's) to estimate, not just 
the next term in the perturbation series, but the entire sum of the series 
(as is frequently done in condensed- matter applications), for both a e and 
a M . We obtain in this way a more accurate theoretical prediction of the 
QED contribution to a M , in particular, which lies outside the errors quoted 
previously. 

The first step is to obtain an accurate value for the fine-structure constant 
a. The two most precise measurements of a areEl 



a' 1 = 137.0359979(32) (7) 



and 
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a' 1 = 137.0359840(50) (8) 

We note that these two values differ by more than 2 standard deviations, but 
nevertheless take the average of eqs (0) and © to obtain 

oi-i = 137.0359939(27) (9) 

The accuracy of this result limits the precision of tests of QED in the case 
of a e , where both theory and experiment are extremely precise. The pertur- 
bation series for a e is@ 

a e = l(-) - 0.328478965(-) 2 + 1.17611(42)(-) 3 - 1.434(138)(-) 4 (10) 

2 7T 7T 7T 7T 

and the error in the theoretical prediction is dominated by the error in a exp . 
The [N/M] Pade Approximant to a series 

S = S + S lX + S 2 x 2 + ... + S N+M x N+M (11) 

is given by 

1 1 J 1 + b lX + ... + b M x M K ' 
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where one chooses the coefficients dj, bj so that 

[N/M] = S + 0(x N+M+1 ) (13) 

One can use such a PA either to predict the next coefficient Sn+m+i or to 
evaluate [N/M] for the relevant value of x (in our case x — - ), and obtain 
an estimate for the sum of the series. Here we do the latter. The PA's are 
known to accelerate the convergence of many series by including the effects 
of hi ghe r (unknown) terms, thus providing a more accurate estimate of the 
seriesli-!J. The PA's also provide reliable estimates of many asymptotic series, 
as is the case in QEE0 and QCdI. 

For our application, we first construct PA's to a e after removing an overall 
multiplicative factor of (-). Our result for the [1/2] PA is 

[1/2] = 1159652169.1(24.0) x 1(T 12 (14) 

and the [2/1] PA agrees very well with the [1/2]: 

[2/1] = 1159652169.0(24.0) x 10~ 12 (15) 

The errors consist of 22.8 from a and 7.4 from the theoretical uncertainty. 
To obtain a e one must add the contribution due to muon diagrams 

Aa e (muon) = 2.8 x 10~ 12 (16) 

the contribution due to r diagrams 

Aa e (tau) = 0.01 x 10~ 12 (17) 

the hadronic contribution 

Aa e (hadron) = 1.6(2) x 10~ 12 (18) 
and the purely weak contribution 

Aa e (weak) = 0.05 x 10" 12 (19) 

for a total of 

Aa e = 4.5(2) x 10~ 12 (20) 
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Thus the theoretical prediction for a e is 



a e = 1159652173.5(24.0) x 10~ 12 (21) 

Comparing eq (^TJ) with eq (|3|), we see that there is beautiful agreement 
between theory and experiment: 

off* - a e = 14.7(24.0) x 10~ 12 (0.61a) (22) 

As noted before, the error in eq is dominated by the error in a. If 
one now assumes that QED is correct, and hence that theory and experiment 
agree, one obtains a new and more accurate value of a: a t h, where the [1/2] 
PA gives 

a th = 137.03599228(86) (23) 

and the [2/1] PA leads to 

at* = 137.03599227(86) (24) 

Comparing eq (|23|) with eq @, one sees that there is beautiful agreement 
with the less-precise experimental value. 

«cxpt - otf = 16(28) x 10- 7 (0.57a) (25) 

corresponding to the good agreement in eq fl22|). We note in passing that 
this provides an a posteriori justification for averaging naively the two most 
accurate measurement of a, and that the difference between the values 

of a^ 1 extracted using the perturbative series and the PA's is just 3 x 10~ 8 . 

We now turn to the anomalous magnetic moment of the muon, a M . As is 
usual, we first consider the difference^'!! 

Ot Ot ex 

a^-a e = 1.09433583(7)(-) 2 + 22.869265(4)(-) 3 + 127.00(41)(-) 4 (26) 

7T 71 71 

In constructing a PA to this series, we must first remove an overall factor 
(^) 2 from the perturbative series. In this way, we obtain the [1/1] PA value 

(a M - a e )[l/l] = 6194839(12) x 10~ 12 (27) 

whereas the value from the series in eq (p6|) is 

-a e = 6194791(12) x 10~ 12 (28) 
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Adding a e from eq (0]), after subtracting Aa e from eq (p0|), we obtain 



ag ED = 1165847008(12)(27) x 10" 12 (29) 

where the first error is due to numerical integrations used in evaluating the 
perturbative series, and the second error is due to a. This should be com- 
pared the value of Kinoshita and MarcianoS 

a|f D = 1165846955(44)(27) x 1(T 12 (30) 

We note that the difference between these two estimates of a« ED is consid- 
erably larger than the error propagated from a. The reason for our smaller 
error is that we have used the new more precise values in ref. [8]. 
If one now adds the hadronic^H and the weak@ contributions 



Aa M (had) = 7011(76) x 10" 11 (31) 

and 

Aa M (weak) = 195(10) x 10~ n (32) 
one obtains the theoretical value 

= 116591907(77) x 10" 11 (33) 

The error is dominated by the error in Aa M (had), and new, more precise 
experiments are underway in Novosibirsk and FrascatS to reduce this error. 
Comparing eq (33) with eq (6), we obtain 



a 



cx p* _ a M = 4(9) x 10~ 9 (0.4(j) (34) 



The error in the difference between theory and experiment is dominated by 
the experimental error in eq (6), which should be reduced by a factor of 20 
in the forthcoming BNL experiment!!. 

In summary, we have used PA to obtain new more precise values for 
the QED values of a e and a M . These PA values, in effect, estimate the 
unknown higher-order contributions, and should be more precise than the 
naive perturbative values used previously. It would be interesting to compare 
our estimates with values obtained in a different way, for example using 
the effective charge approachli^l which agrees very well with PA's in QCD 
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applications^!. Although smaller than some of the other uncertainties, the 
shift we find in a M , in particular, is significantly larger than other theoretical 
uncertainties and the error due to a. 
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